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ABSTRACT 

We report on a measurement of the angular power spectrum of the anisotropy in 
the Cosmic Microwave Background. The anisotropy is measured in 23 different multi- 
pole bands from £ = 54 3°) to i = 404 (^i 0.45°) and in 6 frequency bands from 
26 GHz to 46 GHz over three observing seasons. The measurements are consistent from 
year to year. The frequency spectral index of the fluctuations (measured at low i) is 
consistent with that of the CMB and inconsistent with either dust or Galactic free-free 
emission. Furthermore, the observations of the MSAMl-92 experiment ( |Uheng et al 



1994|) are repeated and confirmed. The angular spectrum shows a distinct rise from 



6Ti = ^i{2e + 1) < |af |2 > /Air = 49t^/iK at £ = 87 to ST^ = 85t^VK at £ = 237. 
These values do not include an overall ±14% (Icr) calibration uncertainty. The analysis 
and possible systematic errors are discussed. 

Subject headings: cosmic microwave background — cosmology: observations 



NRAO, 2015 Ivy Rd., Charlottesville, VA, 22903 



1 



2 



NetterReld et al. 



1. Introduction 

The discovery of the anisotropy in the cosmic mi- 
crowave background (CMB) by the C OBE sateUite 
( Smoot et al. 1992 ; Bennet et al. 1992 ) and its subse- 
quent confirmation ( Ganga et al. 1993 ) at large angu- 
lar scales (greater than 7°) has been followed by pos- 
itive detections of anisotropics in the microwave sky 
at intermediate angular scales by a number of other 
experiments (De Benardis et al. 1994t iDevlin ct al. 



1994; 


Dragovan et al. 1994 


; Gundcrsen ct al. 1995|; 


Gutierrez et al. 1995) and upper limits at smaller 


angular scales ( 


Meyers, Readhead, & Lawrence 1993|; 



Tucker et al. 1993| ). In order to reject foreground con- 



tamination, most of these experiments observe the sky 
over a range of frequencies. Additionally, several have 
successfully repeated their measurements of previous 
years. 

Characterization of the anisotropy at medium an- 
gular scales (2° to 0.2°) can strongly constrain theo- 
ries of structure formation and cosmological param- 
eters, mainly through the angular power spectrum 



of the anisotropy (Bond et al. 1994 


• 


Crittenden & 


Turok 1995 




Kamionkowski et al. 1994 




Jungman et 


il. 1995. For a recent review of the CMB see 


White, 



Scott, fc Silk 1994) 



We present results from an experiment designed 
to measure the angular power spectrum of the CBR 
at medium angular scales. The SK telescope ob- 
serves from the ground in Saskatoon, Saskatchewan, 
Canada. The observing scheme has many internal 
consistency tests which allow checks of the integrity 
of the measurement. A frequency span of 20 GHz 
provides discrimination against foreground contam- 
inants. Measurements made by this experiment are 
consistent from year to year. In addition, we have 
reproduc ed the results of the MSAMl-92 (|Cheng et 
il. 1994| ) experiment. P artial descriptions of the SK 
instrument ar e given in WoUack et al. 1993 (SK93) 
Wollack 199^ (SK93), [WoUack et al. 1994 jPage et&i 



1995 



1994 [Netterfield 1995| (SK94) and [Netterfield et al 
(SK93 to SK94 comparison). The instrument 



and calibration is described in detail in Wollack et al 



1996| . Previous results (SK93 and SK94) are included 
in this analysis. 

2. Instrument and Calibration 

The SK telescope is comprised of a corrugated feed 
horn illuminating a parabolic primary followed by a 
chopping flat with a vertical chopping axis. The tele- 



scope is steerable in azimuth, but fixed in elevation. 
Observations are made in Ka-band (26 GHz to 36 
GHz) and Q-band (36 GHz to 46 GHz). Observa- 
tions made in each band are broken up into 3 fre- 
quency sub-bands and two linear polarizations. To- 
tal power radiometers based on HEMT amplifiers are 
used. The supernova remnant Cas-A is used to cali- 
brate the telescope and determine the beam size and 
pointing. Table |l| lists the beam parameters and un- 
certainties. The beam widths are known to better 
than 2%. The calibration uncertainty is comprised of 
a 13% contribution due to uncertainty in the absolute 
calibration of Cas-A, which is constant between years 
and radiometers, and a 3% to 5% contribution due to 
measurement uncertainty. 

The telescope pointing is determined from the po- 
sition of Cas-A. In 1995, the telescope was pointed 
0.05° differently than was intended, which was in- 
cluded in the analysis. In 1994, the discrepancy was 
less than 1/10 of a beam, and was ignored. In addi- 
tion, there was a 0.03° (Itr) jitter in the beam position 
from the chopping plate and base pointing accuracy. 
When convolved with the beam, this can be approx- 
imated by adding this uncertainty in quadrature to 
the nominal beam width. For 1995 this is results in a 
1% widening of the effective azimuth beam. 



3. Observing Strategy: Synthesized Beams 

The observing strategy is designed to offer a vari- 
ety of internal systematic checks, minimize contam- 
ination by atmospheric temperature gradients, and 
to simultaneously probe a variety of angular scales. 
These goals are achieved by sweeping the beam in az- 
imuth on the sky by many beam widths with a large 
chopping flat, and then synthesizing effective antenna 
patterns in software. For simplicity in the discussion 
which follows, telescope and observing parameters for 
the 1995 season are used. Parameters for all years are 
listed in Table 0. 

As the beam is swept on the sky, the radiometer 
is sampled 168 times per complete sweep. In analysis 
the 168 samples from each sweep, T^, are multiplied 
by a weighting vector, Wi, to give ATg^eep = J2 "^iTi- 
By choice of weighting vector the relative weighting 
of each spatial point in the sweep can be set allowing 
the synthesis of arbitrary effective antenna sensitivity 
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Table 1 
Observing Parameters 







Ka93 


Ka94 


Q94 


Q95 Cap 


Q95 RING 


Beam 


X FWHM 


1.42 ± 0.005° 


1.42 ± 0.005° 


1.004 ± 0.005° 


0.461 ± 0.003° 


0.461 ± 0.003° 




y FWHM 


1.42 ±0.02° 


1.42 ±0.02° 


1.08 ±0.02° 


0.513 ±0.004° 


0.513 ±0.004° 


Chopper 


Pattern 


Sine 


Linear 


Linear 


Sine 


Sine 




Amplitude on Sky 


4.90° 


7.00° 


7.35° 


7.96° 


3.36° 




Rate 


3.906 Hz 


3.906 Hz 


3.906 Hz 


2.976 Hz 


2.976 Hz 




Samples / Sweep 


16 


64 


64 


168 


168 


Pointing 


Azimuth (W:E) 


-8.2° : 7.8° 


-7.15° : 7.25° 


-7.15° : 7.25° 


-7.392° : 7.288° 


0.14° 




Jitter 




0.02° 


0.02° 


0.03° 


0.03° 




Elevation 


52.2° 


52.2° 


52.2° 


52.24° 


60.2° 


Timing 


Time Per Base Move 


16.4 s 


20 s 


20 s 


40.3 s 


2 weeks 




File Length 


17 min 


15 min 


15 min 


20 min 


20 min 



patterns, H{x,y), given by 

iJ(x) = /^m,G(x-X0\ (1) 

\ ^ / bin 

where is the position on the sky of the center of the 
main beam corresponding to weighting vector element 
Wi and 

G{x, y) = exp f-^ - ^) (2) 

is the main beam pattern of the telescope pointed at 
(0,0). The azimuth dimension is denoted by x and 
the elevation dimension by y. The beam width of the 
telescope is given by cr = FWHM/^S ln(2). The ef- 
fect of the rotation of the Earth during an integration 
is included in i?(x), as denoted by < . . . >un- 

For example, if the weighting vector wi is selected 
so that data which are acquired when the chopping 
plate is centered are given a positive weight, and data 
which are acquired when the chopping plate is offset 
to the left or right are given negative weight, then 
a three lobed synthesized beam on the sky is pro- 
duced (hereafter referred to as the 3pt synthesized 
beam). The weighting vectors are normalized so that 
J |i/(x)|(ix = 2 for a single sweep, not including the 
rotation of the Earth. The mean of the weighting 
vector is zero in order to eliminate sensitivity to fluc- 
tuations in the mean temperature of the sky. Since 



the weighting vectors for different effective antenna 
patterns are made to be nearly orthogonal all of the 
effective antenna patterns may be acquired from the 
same data. Synthesized antenna patterns up to 19pt 
(« 60 Hz) are produced. The number of lobes is lim- 
ited by the beam size to throw ratio, and not by the 
detectors. 

The synthesis of the 9pt beams is demonstrated in 
Figure |l|. Figure § displays some of the synthesized 
beams. In total, including the three years, there are 
37. The 3pt through 9pt antenna patterns are sim- 
ilar for Ka94, Q94 and Q95, allowing these data to 
be compared. Since the sweep amplitude and point- 
ing for SK93 was substantially different than for the 
other years, the Ka93 3pt beam is also substantially 
different. Comparing SK93 with SK94 requires the 
use of a specially optimized weighting vector (Ka93 
overlap). 

Data are acquired in two modes: In the first mode 
[Cap Data), which was used all three years, the chop- 
ping plate axis is vertical which places the beam at 
a constant elevation of 52.2°, the elevation of the 
North Celestial Pole (NCP) in Saskatoon. The beam 
is swept sinusoidally 12.8° in azimuth (7.35° peak to 
peak on the sky) at 2.976 Hz. The center of the sweep 
is alternated every 40.3 seconds between 7.32° East 
of the NCP and 7.28° West of the NCP. As the earth 
rotates throughout the day, the celestial polar cap is 
covered, from the NCP to 82° declination. Data ac- 
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Q95 9pt 
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Fig. 1. — Beam Synthesis. Shown is the synthesis of 
the 9pt beams for Q94 and Q95 using equation (Q). 
H{x) is shown along the path of the swept beam. The 
effect of binning in RA is neglected here for clarity. 
In 1994, an approximately linear sweep pattern, Xj, 
was used as shown. This allowed the use of sinu- 
soidal weighting vectors, Wi, to produce synthesized 
beams with equally spaced lobes. The use of a sinu- 
soidal sweep rendered this impossible in 1995. The 
Q95 weighting vectors were generated by optimizing 
the sensitivity of the synthesized beam to a specified 
band of angular scales. Spatially similar synthesized 
beams are produced. 
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Fig. 2. — Synthesized Beams. The contour lines are 
every 3dB. Dashed lines denote negative weighting. 
The NCP is located at (0,0). Straight lines from the 
NCP are lines of constant hour angle. The bottom 
right plot is of the Ring 4pt synthesized beam, which 
is centered 8° above the NCP. The remaining plots 
are of Cap data acquired in the East. The curvature 
is due to the fact that the beam is swept in azimuth, 
and not declination. The lobes farthest from the NCP 
are wider in RA than those closest because the data 
are averaged into 24 (3pt to 5pt) or 48 RA bins. 
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quired in the East are repc;ated « 12 hours later in the 
West, giving a very powerful systematic test. These 
data are multiplied by weighting vectors which pro- 
duce a family of synthesized beams on the sky from 
3pt to 19pt, and are then integrated into 24 (for 3pt 
to 5pt) or 48 (for 6pt to 19pt) bins in Right Ascension 
(RA). 

In the second mode {Ring Data), which was only 
used in 1995, the chopping plate is tilted back 4.0° 
which raises the beam elevation 8.0° to 60.2°. The 
telescope is pointed so that the center of the sweep is 
North. The beam is swept sinusoidally with a peak 
to peak amplitude of 3.36°on the sky. In the small 
angle approximation, the beam is swept in RA at a 
constant declination. As the Earth rotates an entire 
ring at 82° declination is sampled. 

Synthesized beams arc generated in the same way 
as with the Cap data, except that a new weighting 
vector must be generated every sweep to keep the 
synthesized beam fixed in RA during its integration. 
The synthesized beams are smaller than the ampli- 
tude of the sweep which allows positions on the sky 
to be tracked in software. Three point to 6pt syn- 
thesized beams are produced, as well as the effective 
3pt beam pattern of the MSAMl-92 experiment. The 
Ring data lacks the East to West comparison test of 
the Cap data. However, synthesis of the MSAM beam 
makes it possible to compare directly with a very dif- 
ferent experiment. 
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4. Data Selection and Reduction 

Because the SK telescope observes from the ground 
in Saskatoon, an important part of the analysis is 
the selection of data where contamination from atmo- 
spheric spatial temperature fluctuations do not domi- 
nate. The first cut is made by evaluating, for each 20 
minute file, the mean deviation of 8 second averages 
made with a beam synthesized using a 2pt weighting 
vector which is sensitive to the horizontal component 
of spatial temperature gradients. This number is large 
when the spatial gradient is changing and small when 
it is stable and has been found to be a good indicator 
of atmospheric contamination. The mean deviation 
is used because it is less sensitive than the standard 
deviation to spikes in the data, which are more likely 
to be due to birds or airplanes than to atmospheric 
noise. The mean deviation is related to the standard 
deviation by A^, = ^^<^x for normal distributions. 

A more general expression for the cut levels is given 



Fig. 3. — The sky noise in Saskatoon. The value (" is 
a measure of the stability of the horizontal component 
of spatial temperature gradients (see text). The upper 
dashed line in each plot denotes the nominal atmospheric 
cut. The lower dashed line indicates the system noise 
limit, Cn = KNET/Oeff. The noise floor for the 1995 
Ring data is higher than for the Cap data because O^f / is 
smaller. Data with ^ > 10 mK- s^'''^/deg have been lim- 
ited to C = 10 mK- s^^^/ deg in this plot. The maximum 
value of C recorded was 1360 mK- s^/^/deg. 
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by 

C - <jV^/9eff (3) 

where a is the standard deviation of the 2pt data, r 
is the integration time, and 

Oeff — J dxH{0x,9y)d9xd0y (4) 

is the effective sweep amphtude in degrees. The units 
of ^ are mK ■ s^^^ / Aeg. Neglecting the radiometer 
contribution, this measurement of the stabihty of the 
atmosphere may be compared from experiment to ex- 
periment. The effective sky noise, C,n, due to the 
receiver NET is given by C„ = KNET/6'e//, where 

~ v^^^Ci^^' ^"^^ ^ number of samples per 

sweep. 

Figure ^ shows C as a function of time for each of 
the observing seasons. Horizontal dashed lines are the 
nominal cut levels (the upper line) and the receiver 
noise level, C„ (the lower line), for each of the observ- 
ing runs. Each point represents the result for one file 
of data (see Table p. Because the throw for the Ring 
data was significantly smaller than for the Cap data, 
Qn is significantly higher. In other words, the Ring 
observations are significantly less sensitive to the at- 
mospheric gradients than the Cap observations. 

It is furthermore required that for a file to be ac- 
cepted, the files before and after must also pass the 
atmospheric cuts. Thus, the sky must be stable for 
three files (45 minutes in 1994 and 60 minutes for 
1995) for any data to be accepted at all. Of the data 
acquired, 27% passes the nominal cut. Table || sum- 
marizes the amount of data acquired in each observing 
mode at different cut levels. The cut level listed first 
for each experiment is the one used for the definitive 
data set. The lower cut data for the Ring is com- 
prised of the two long receiver noise limited stretches 
around day 59 and day 63 of 1995 rather than be- 
ing based on specific sky noise level. The analysis is 
performed for each of the listed cut levels (see Ta- 
ble |ll|). Decreasing the cut level from the definitive 
values changes the spectrum by approximately la for 
the SK95 data. Because this lower cut level is so close 
to the receiver noise, good data are cut at random, 
reducing the length of contiguous stretches of data. 
This increases sensitivity to long time scale noise and 
reduces the robustness of error bar calculation. In- 
creasing the cut level has an insignificant effect on 
the results. 

The data are blanked for base moves (e.g. 4 sec- 
onds every 40 seconds for SK95) and for spikes which 



may be due to birds, planes or occasional glitches in 
the data system. An entire sweep of data is blanked 
if any of the samples in it (64 for SK94, and 168 for 
SK95) exceed the mean by more than some limit (3.5cr 
for SK94 and 3.85cr for SK95). The amount of data 
removed in this step is consistent with the data being 
normally distributed (15% for SK95). 

The data are then multiplied by the weighting vec- 
tor and binned according to RA. A weighted mean of 
the frequency and polarization channels is also made 
before application of the weighting vector for analyses 
where frequency and polarization information are not 
required. For the Ring data, a new weighting vector 
is generated every sweep so that the antenna pattern 
in each bin is fixed in RA, while for the Cap data, the 
weighting vectors are fixed. The Cap data acquired 
in the East are binned separately from the data ac- 
quired in the West. The bin numbering for both East 
and West data is based on the RA of the center of the 
sweep when the telescope is pointed in the East. 

For the Cap data, the error bars are generated in 
two stages. First, after beam synthesis, the data are 
averaged into bins 15 sidereal minutes wide. Weights 
for each bin are calculated based on the distribution 
of 18 second averages. If any 18 second average is 
more than 3.5cr away from the mean of the bin it is 
cut. These 15 min averages are then binned into 24 
or 48 bins around the sky with weights propagated 
from the weights of the 15 minute averages. These 
weights are used to give the relative weighting of each 
RA bin. The overall weighting is generated from the 
distribution of the 15 minute bins after removal of cor- 
responding mean signal from the 24 or 48 RA bins. 
For data weighted to synthesize many lobed beams 
(6pt or more) the error bars based on 18 second av- 
erages agree to 10% with those based on 15 minute 
averages, but for the 3pt or 4pt data (which is more 
sensitive to atmospheric noise) the distribution of 15 
minute averages predicts error bars as much as 50% 
larger than would be generated from the distribution 
of the 18 second averages (see Table ||). For the Ring 
data, with its substantially smaller sweep which yields 
much less sensitivity to the atmosphere, the error bars 
are generated directly from the distribution of the 18 
second averages. 

5. Offsets 

The design of the SK telescope is intended to min- 
imize the instrumental offset. However, the emis- 
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sion from the large chopping plate is a function of 
the plate orientation, which theoretically produces a 
~ —400 offset for the horizontally polarized chan- 
nel and a « 200 fiK offset in the vertically polarized 
channel in the 3pt data. The measured 3pt offsets are 
approximately this magnitude for SK93 and SK94. 
In SK95, however, an additional w — 1 mK offset was 
present in both polarizations. The source of this ad- 
ditional offset is not known. See WoUack ct al. 1996 
for possible explanations. 

The offsets are dealt with differently between the 
Cap and Ring data analysis. In the Cap analysis, the 
weighting Wi is constant for a given synthesized beam 
and chopping plate orientation. For this reason the 
offsets are ignored until the data have been multi- 
plied by the weighting vector and binned in RA. At 
that point, the mean value of the 24 or 48 RA bins 
is removed. The error bars are then multiplied by 
y'24/23 or \/48/47 to compensate for the removed 
degree of freedom. The levels of the Cap offsets for 
each synthesized beam are listed in Table ||. 

For the Ring data, however, a new weighting vector 
must be calculated for each sweep, which causes the 
offset to change as the synthesized beam is tracked 
over the sky. To deal with this, all of the accepted 
data are coadded synchronously with the sweep to 
form a 168 sample mean sweep. This mean sweep, 
which is essentially the instrumental offset as a func- 
tion of chopping plate position, is subtracted from 
the data for each sweep prior to multiplication by the 
weighting vector. At this point then, the offset for 
any synthesized Ring beam is zero. 

As evidenced by the larger in scatter of the 15 
minute averages compared to that predicted by the 
distribution of the 18 second averages in the Cap data 
(see Table ^), there may be some drift to this offset 
on longer time scales. Since at least some of offset 
is due to emission, the excess long time scale noise 
could be expected to correlate with the temperatures 
of the optical components. The correlations are not 
significant la) however, and regressing out the op- 
tics temperatures has no effect on the long time scale 
noise. Since the gain of the amplifiers is a function 
of temperature the long time scale noise might be ex- 
pected to vary with receiver temperature. But, as 
with the optics temperatures, there is no significant 
correlation with receiver temperature, and regressing 
it out has no effect. In the final analysis, nothing is 
regressed out. 

To further investigate the effect of the excess long 



time scale noise in the Cap data, polynomials in time 
(linear to 9th degree) were fit out of the data. The 
maximum slope from fitting a line in time to the data 
is 7 fiK per day for the SK95 3pt data. The typical 
slope is less than 1 fiK per day and consistent with 
zero. As with regressing out the auxiliary tempera- 
tures, removing lines or polynomials has no significant 
effect on the long time scale noise or on the final an- 
gular power spectrum. In the final analysis, only the 
constant offset is removed. 

To test the possibility that the long time scale noise 
is due to instrumental drifts, data taken in the West 
can be subtracted from data taken in the East 40 sec- 
onds earlier. Doing this decreases the differences be- 
tween the short and long time-scale based error bars. 
This subtraction, however, has no significant effect 
on the spectrum, beyond increasing the error bars: 
{E - W) / V2 will have the same rms, as i? or ly by 
themselves, but now there is only one data set to do 
statistics on, rather than two. It is also possible that 
the excess long time scale noise is due to atmospheric 
contamination: this was the interpretation taken by 
the OVRO Ring experiment ( Readhead et al. 1989[ ) 
for a similar effect. 

6. Internal Consistency Checks 

The observing scheme used in the SK experiment 
has a number of internal consistency checks. For the 
Cap data, the rotation of the Earth causes the same 
sky as was observed with the telescope pointed to the 
East of the NCP to be observed « 12 hours later 
when the telescope is pointed to the West. This test 
{E — W) gives assurance that the signal which has 
been observed is truly on the sky: contamination from 
the atmosphere or ground pickup could be expected to 
repeat on a 24 solar hour time scale, not a 24 sidereal 
hour time scale with the signal in the East lagged 12 
sidereal hours from the signal in the West. 

The data set may also be divided into first half 
(HI) and second half {H2), and compared. Since 
there is approximately one month between the cen- 
troids of the first and second half, the success of this 
test gives assurances that the observed signal repeats 
on sidereal time, and not solar time: if the signal re- 
peated on solar time, the signal between halves would 
be lagged by 4 out of 48 bins. This places hmits on 
contamination by the sun and any Radio Frequency 
Interference (RFI) that recurs daily. 

The data from the two polarization channels {A 
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Table 4 

Internal Consistency Tests I: 3pt to 9pt Cap 





Ka93 Overlap 


3pt 


4pt 


5pt 


6pt 


7pt 


8pt 


9pt 


Ka93 Sum 


3.76 
















E-W 


1.17 
















A-B 


1.0 
















Ka94 Sum 


5.20 


4.59 


2.16 


2.92 


1.47 


1.40 


1.13 


1.03 


A-B 


0.79 


0.55 


0.67 


0.70 


0.85 


0.94 


0.88 


1.18 


E-W 


0.76 


0.88 


0.83 


0.70 


1.17 


1.06 


1.39 


0.59 


Quadrature 




1.2 


0.7 


1.0 


1.2 


1.3 






Q94 Sum 


3.21 


2.15 


1.97 


1.14 


1.06 


1.27 


1.19 


1.53 


E-W 


0.91 


0.67 


1.06 


1.09 


0.68 


1.08 


1.05 


0.84 


Quadrature 




0.7 


1.1 


1.0 


1.1 


1.4 






Q95 Sum 




2.64 


1.98 


3.27 


2.49 


2.49 


2.69 


2.37 


A-B 




0.91 


1.05 


1.09 


0.74 


0.73 


1.1 


1.0 


E-W 




1.13 


0.99 


1.22 


1.12 


1.60 


1.09 


1.16 


H1-H2 




2.34 


1.60 


1.05 


1.13 


1.09 


1.06 


0.88 


Quadrature 




0.93 


1.23 


1.37 


1.49 


0.90 


1.21 


0.99 


Ka93+Ka94 


6.12 
















Ka93-Ka94 


1.18 
















Ka94+Q94 


7.17 


6.01 


2.91 


3.35 


1.66 


1.52 


1.50 


1.58 


Ka94-Q94 


1.25 


0.73 


1.22 


0.71 


0.87 


1.15 


0.82 


1.03 


Ka94+Q95 




6.01 


2.94 


5.21 


2.53 


2.50 


2.66 


2.56 


Ka94-Q95 




1.22 


1.20 


0.98 


1.43 


1.29 


1.17 


0.85 


Q94+Q95 




4.10 


2.85 


2.92 


2.07 


2.69 


3.10 


2.94 


Q94-Q95 




0.69 


1.09 


1.49 


1.48 


1.06 


0.78 


1.02 



Note. — Reduced for a variety of internal symmetry tests for 3pt to 9pt data. The Sum entries indicate the 
reduced for the total signal on the sky (with A and B polarizations averaged and data acquired in the East and 
in the West averaged). The E — W entries indicate the residual signal if data acquired in the West is subtracted 
from data acquired in the East 12 hours earlier. The A — B entries represent the residual if the two polarizations 
arc difFcrcnccd. The Quadrature entries give the amplitude of the signal if the weighting vector is modified to have 
no sensitivity on the sky. The K^^QS Overlap, 3pt, Apt and 5pt entries have 23 degrees of freedom. The rest have 
47 degrees of freedom. The Sum entries have consistently larger reduced entries than the the control entries, 
indicating the presence of a true sky signal. 
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Table 5 

Internal Consistency Tests II: IOpt to 19pt Cap 





lOpt 


llpt 


12pt 


13pt 


14pt 


15pt 


16pt 


17pt 


18pt 


19pt 


Q95 Sum 


1.95 


2.12 


2.29 


2.30 


1.71 


1.94 


1.62 


1.70 


1.51 


1.18 


A-B 


1.01 


1.36 


0.78 


0.85 


1.36 


0.96 


0.88 


1.04 


1.17 


0.97 


E-W 


0.85 


0.81 


1.29 


1.42 


0.90 


0.85 


1.21 


1.21 


1.01 


1.36 


H1-H2 


0.91 


1.01 


1.29 


0.91 


1.53 


1.62 


0.73 


1.29 


1.19 


1.57 


Quadrature 


0.97 


0.90 


1.22 


1.05 


1.30 


1.31 


0.81 


1.14 


1.09 


1.26 



Note. — Reduced for a variety of internal symmetry tests for lOpt to 19pt SK95 Cap data. All entries have 
47 degrees of freedom. 



Table 6 

Internal Consistency Tests III: Ring 



.'-Spl 4i)t r)pt Opt MSAM 0\('rlap 



Q95Sum 1.86 1.83 1.34 1.01 1.87 

A-B 1.16 1.31 1.40 1.28 1.55 

H1-H2 1.19 1.07 1.24 1.2 1.05 

Q95+MSAM1-92 ••• ••• ••• ••• 3.43 

Q95-MSAM1-92 ••• ••• ••• ••• 1.05 



Note. — Reduced ^oi a variety of internal symmetry tests for SK95 Ring data. The 3pt to 6pt entries have 
95 degrees of freedom. The MSAM Overlap entries have 80 degrees of freedom. 
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is vertical and B is horizontal) may also be com- 
pared. This test is useful for detecting contaminants 
which effect the two receiver chains differently, such 
as data system pickup or polarized diffraction off of 
the ground screen. 

Since the same region of sky has been observed in 
1994 as 1995, the 3pt to 9pt Cap data can also be com- 
pared between years. A comparison of the SK93 and 
SK94 data using a specially prepared weighting vector 
for SK94 is presented in Netterfield et al. (1995) and 
repeated here. Since many elements of the telescope, 
radiometer, and ground shield were changed between 
years, this test can give assurances that the data have 
not been contaminated by side lobes, vibrations or at- 
mospheric gradients. 

Finally, weighting vectors may be applied which 
have the same frequency as sky-sensitive weighting 
vectors, but with the phase adjusted to have no sen- 
sitivity to the sky (quadrature phase). This test is 
sensitive to synchronous vibrational or electrical sig- 
nals, which would not necessarily have the same phase 
as a true sky signal. 

The results of this set of internal consistency tests 
are given in Tables ^||. The expected reduced 
for the difference tests should be distributed around 
1.0 with a standard deviation of l/-\/2z/, where v is 
the number of degrees of freedom, assuming that the 
sky signals are 100% correlated and the error bars 
are completely uncorrelated. This is the case for the 
HI — H2 tests, but not for the E —W tests, the 
A — B tests, or the year to year tests. The small 
beam, large throw in azimuth (rather than declina- 
tion) and slight asymmetry in the East to West point- 
ing of the SK95 experiment means that the E — W 
lagged correlation is not 100%. For the 12pt data, for 
instance, the theoretical E — W correlation coefficient 

(see Section 0) is pf ^ e 

For the 3pt data, ^ = 0.90. In this case, where 
the sky signal is not 100% correlated, we findPI, 



r<EW 



/^Cl^Cf ^ = 0.75. 



2p+(l-p)x|;+w/^ 



(5) 



^For two correlated zero mean normal deviates {x} and {y} with 
(T^ = cr^ , we write yi = ai + pxi where {a} is uncorrelated with 
{x} and p is the correlation coefficient between {x} and {y}. 
One finds that c^+j, = (1 + p)"-'^ and = (1 — p)crx- After 

identifying {x} with the East data and {y} with the West data, 
and noting that x^/*^ = ^~^'^x±y/'^'^ (where a is the instrument 
noise), we find equation (a). 



For the 12pt data, this results in Xe-w/'^ ~ 1-19±0.1 
which is consistent with the value of 1.29 which was 
measured from the data. Similarly, atmospheric noise 
correlates the error bars from the A and B polar- 
izations in the large angular scale data. In this 
case we expect x'a-b/'^ = (1 ~ Pd^)j where = 
/ \J C^^C^^ is the noise correlation coefficient. 

While the majority of the entries in Tables ^ to 
^ are consistent with the the signal being fixed on 
the sky, there are some entries which indicate trou- 
ble. The worst entry is for Q95 3pt HI - H2. The 
value of the reduced is 2.34, which is inconsistent 
with the expected value of 1.0. The HI — H2 resid- 
ual is dominated by long drifts (w 12 hr), which is 
consistent with the earlier observation that the 3pt 
data are especially effected by excess long time scale 
noise. However, the E — W tests, and more signifi- 
cantly, the year to year tests work well for this data 
set, indicating that if the entire season is used the ef- 
fect of the drifts is averaged out. A similar effect at a 
lower level is seen in the Q95 4pt data. Without the 
ability to compare with the other seasons we would be 
much less sanguine regarding the quality of the Q95 
3pt and 4pt data. 

Other notably inconsistent tests are the 14pt and 
15pt HI - H2 tests. For both of these the E-W 
tests are consistent with a true sky signal. However, 
without a comparison between years, we have no ex- 
planation for this. Despite these problems there is 
strong evidence that the measured signal is on the 
sky. 

7. Likelihood Analysis 

In the analysis of these data, two questions are 
asked: 1) is the spectral index of the fluctuations con- 
sistent with that of the CMB, and 2) what is the angu- 
lar power spectrum of the fluctuations. Comparison 
of the data with theoretical predictions may be done 
by comparison with the angular power spectrum. The 
determination of parameters (such as the amplitude 
of the fluctuations at various angular scales, and their 
spectral indices) are made using maximum likelihood 
tests. 

The Likelihood (or the probability density of the 
data, assuming a theory) is defined as 



^- (27r)w/2|M|i/2 
where t is a data set and M = Cjj 



(6) 

Ct is a cor- 
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responding covariance matrix. Cij is the data co- 
variance matrix, and Ct is the pixel to pixel theory 
covariance matrix which is a function of the model 
parameters. The model parameters are adjusted to 
maximize L. For the Saskatoon analysis, the elements 
of the data vector ti correspond to data taken look- 
ing at a specified RA bin on the sky, with a specified 
synthesized antenna pattern, in a given frequency and 
polarization channel. 

Confidence intervals on a parameter T arc found 



from the distribution of L{T). The best value for 
the parameter is determined by finding Tmax where 
L{Tmax) is maximized. For positive definite quanti- 
ties (such as the amplitude of the fluctuations) a non- 
zero detection is claimed when L{Q) / L{Tmax) > 0.15. 
In this case the error bars are found by finding 
such that 



/;-L(r)dT 



= 0.1587 



lo LiT)dT 
for the lower limit, and finding so 



L{T)dT 
j^L{T)dT 



= 0.8413 



(7) 



(8) 



for the upper limit. When L{Q) / L{Tjnax) < 0.15, the 
95% upper limit is found by finding T' such that 



Io'HT)dT 
j^L{T)dT 



0.95. 



(9) 



For non positive definite quantities (such as spectral 
indices) the technique is similar, except that the lower 
limits of the above integrals are —oo rather than 0. 

The data covariance matrix, Co, describes the 
signal due to instrumental noise in the data. The 
diagonal matrix elements are simply the variance 
of the data derived from the distribution of the 15 
minute averages as described in Section 4. The off- 
diagonal elements, summarized in Table |7[ describe 
other instrumental correlations. For instance, the 
noise between frequency channels in the same ra- 
diometer chain is correlated by gain fluctuations in 
the HEMT amplifiers ( |Jarosik 199^ ). Similarly, at- 
mospheric noise introduces correlations between si- 
multaneously acquired channels. Since the weighting 
vectors for SK95 are not strictly orthogonal in time, 
there will be noise correlations between, for example, 
the 13pt data and the 15pt data. The noise from data 
acquired in the East is also correlated with the noise 
from data acquired in the West 20 seconds earlier. 



This is related to the excess long time-scale noise dis- 
cussed previously, and may be due to the atmosphere, 
or drifting instrumental offsets. The noise correlation 
coefficients are generated from the distribution of 18 
second synthesized beam averages for all simultane- 
ously acquired data. All known correlations are in- 
cluded in the analysis. 

The theory covariance matrix, Cy, describes the 



signal due to a hypothesized sky signal ( Bond 1995| 
and references therein; Peebles 1994; White & Scott 



1994). With the theoretical temperature fluctuations 



on the sky expressed in spherical harmonics as 



7^(x) = Ea^' 



(10) 



the theory covariance matrix can be expressed as 

=^E(2^+1)^^W^^'' (11) 
where = (|a™P) ^'^'^ 

Wl' = Jditi y"dx2i?.(xi)i7,(x2)PHxi •X2) (12) 

is the window function. i?(x) is the effective antenna 
pattern from equation (|l|), and Pi{xi ■ X2) are the 
Legendre polynomials. The window functions associ- 
ated with the Saskatoon experiment are presented in 
Figure ||. 

Rather than expressing the sky fluctuations in 
terms of Ci, we use 



bit = \l Q. 

47r 



(13) 



In terms of 5Ti the theory covariance matrix becomes 

Ml' 



(14) 



So ST^ is the variance per logarithmic interval of the 
angular spectrum of the CMB. 

8. The Window Function 

The angular scale to which a data set is sensitive 
{£e) is determined from the diagonal components of 
the window function, W^" as follows. The expected 
rms amplitude of the fluctuations. A*, predicted by 
a theory for data acquired with a specified beam is 
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Table 7 
The Correlations 



Between Source Magnitude ( 



Frequency Channels HEMT correlations, Atmosphere 0.2 - 0.7 

Polarization Channels Atmosphere 0.0 - 0.3 

East /West Atmosphere, Offset drift 0.0 - 0.2 

Q95 npt and (n + 2)pt non-orthogonal weighting vectors 0.0 - 0.2 



Note. — A summary of the noise correlations included in the analysis and their sources. 



K,/9eff Cut Level C hours 
(mK-sVydeg) 



Ka93 


0.60 


4.5 


130 






3.0 


79 


Ka94 


0.44 


2.5 


140 






1.7 


107 


Q94 


0.42 


2.5 


158 






1.7 


98 


Q95 Cap 


0.36 


2.5 


174 






3.5 


242 






1.9 


114 


Q95 Ring 


0.86 


6.8 


121 






9.7 


132 






By Hand 


66 



Table 2: Amount of data at different cut levels. The 

first cut level for each experiment listed is the nominal 
cut used. The others are used for consistency checks. 
K/6eff is used to convert ( to receiver NET (see text). 




Fig. 4. — The window functions for each synthesized 
beam for SK93-SK95. The many-lobed Cap window 
functions are narrower than the Ring windows, pro- 
viding finer angular resolution. 



The Angular Power Spectrum of the CMB 



13 







A Offset 


B Offset 


Ka93 3pt 


1.2 


-60 ± 3 /iK 


-360±3/iK 


Ka94 3pt 


1.1 


117±6 


-550 ±7 


4pt 


1.3 


10 ±5 


-50 ±6 


5pt 


1.0 


54 ±6 


225 ±6 


6pt 


1.0 


1±7 


5±7 


7pt 


1.0 


7±9 


62 ±9 


8pt 


1.0 


0± 13 


34± 13 


9pt 


1.0 


80± 19 


97± 18 


Q94 3pt 


1.2 




-695 ±8 


4pt 


1.2 




63 ±7 


5pt 


1.0 




281 ±8 


6pt 


1.0 




15 ±8 


7pt 


1.1 




79 ±9 


8pt 


1.0 




1 ± 11 


9pt 


1.0 




51 ±13 


Q95 3pt 


1.5 


-616 ± 13 


-2182 ±9 


4pt 


1.4 


4± 11 


31 ±7 


5pt 


1.2 


89± 10 


436 ±6 


6pt 


1.2 


63 ±9 


107 ±6 


7pt 


1.1 


21 ±9 


145 ±6 


8pt 


1.2 


-56 ±9 


-83 ±6 


9pt 


1.1 


21 ±9 


53 ±6 


lOpt 


1.2 


-101 ±10 


-103 ±6 


llpt 


1.1 


25± 10 


35 ±6 


12pt 


1.1 


70± 10 


75 ±6 


13pt 


1.1 


-26± 11 


7±7 


14pt 


1.0 


-34± 12 


-41 ±7 


15pt 


1.1 


-23± 12 


-31 ±8 


16pt 


1.1 


-12± 13 


-17±8 


17pt 


1.0 


15± 14 


12 ±9 


18pt 


1.1 


-41 ±15 


-49 ±9 


19pt 


1.0 


-65± 17 


-31 ±10 



Table 3: Offsets and long time scale noise. The ratio 
of error bars based on the distribution of 18s aver- 
ages to error bars based on the distribution of 15 min 
averages is crismm/ciss- The 15 minute averages are 
used. A Offset is the offset in the vertically polarized 
channel. B Offset is the offset in the horizontally po- 
larized channel. For Q94, the A receiver chain was 
damaged in shipping and was not used. 



given by the theoretical correlation matrix (eq. [|14|) 
as 



(15) 



or the sum of the theoretical angular spectrum, weighted 
by W" /i. Under the assumption of a flat spectrum 
in 5Tg,, this may be simplified to 



A' = 5T, 



where 



(16) 



(17) 



This provides a convenient conversion between the 
rms amplitude of the fluctuations, A, and the flat 
band power estimation, 5Ti^. The effective ^ for a 
single window function is the centroid of Wl^ /I, or 



(18) 



The band power estimate, 5TI , over a window 
function, W" , can be generated directly from the an- 
gular spectrum, STg, by 



sn 



(19) 



When combining the results of several measurements 
with differing window functions it is appropriate to 
find the weighted mean of ST^i , since these are pro- 
portional to the variances, which are normally dis- 
tributed far large degrees of freedom. Thus, 



6T, = 



(20) 



Inserting equation ([T^) into equation (|2^) we find 



ST, 



where we have introduced the definition 



(21) 



(22) 



Wi is the effective window function which results 
from combining band power measurements with dif- 
ferent window functions. The quantity (J^t^ is the 

measurement uncertainty on ST'^i . Note that / [W = 
1 with this definition. 
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9. The Theoretical Power Spectrum 



associated with this is 



If the angular spectrum varies slowly over the an- 
gular scales to which the experiment is sensitive it 
is reasonable to expand STg as a power law in both 
angular scale, £, and observing frequency, ly, as 



(23) 



where STi^ is the amplitude of the fluctuations at an- 
gular scale £e and frequency i^o- The angular spec- 
tral index is m and the frequency spectral index is f3. 
With £e chosen using equation (^8|) and (3 chosen as 
the center frequency of the observations (i.e. ,36 GHz) 
the three parameters are roughly orthogonal for m 
and P near 0. Therefore, neither the amplitude of the 
fluctuations, ST{, nor the frequency spectral index, /3 
will be a function of choice of angular spectral index, 
m. 

Determination of the frequency spectral index, /3 
is useful for discriminating the source of the fluctu- 
ations. For the CMB /3 = 0, while in the Rayleigh- 
Jeans region for Galactic free- free emission (3 — —2.1, 
and for Galactic dust (3 = 1.7. 

By combining data from multiple window functions 
the angular spectral index, m, may be determined. 
Since the limits on this parameter are often uninter- 
estingly large, it is often simply set to some fixed 
value. Choosing to = assumes constant STi over 
the region where the data are sensitive. The deter- 
mination of 5Ti^ for TO = is known as a flat band 
power estimate of the angular spectrum. The choice 
of m = 1 reproduces the "delta function" correlation 



3t al. 1995 



function used in WoUack et al. 1993 and Netterfield 



The procedure for finding the angular spectrum, 
5Ti, is to group data sensitive to a particular range 
of angular scales (e.g., the 5pt data are sensitive to 
angular scales around £ = 100) and then use the max- 
imum likelihood test to find the best amplitude, ST^^, 
and where possible the spectral indexes, P and m, at 
each angular scale. Where there is insufficient lever- 
age to determine /3 and m, they are fixed at 0. 

Of mainly historical interest is the Gaussian auto- 
correlation function. 



Ci?(ei2) = C(0) / dxi / dx2H,(xi)//j(x2)oxp 



2^2 



(24) 

where 9i2 — acos(xi • X2)- The angular spectrum 



ST, 



'c(o,m»p 



2£! 



(25) 



with £o = 1/Oc- The coherence angle of the theory is 
9c, and is often varied to maximize overlap with the 
experiment. The free parameter which is varied in a 
maximum likelihood test is C(0). This spectrum has 
no particular theoretical motivation. However, since 
it changes only slowly over the window function of 
most experiments, C(0) can be converted to a band 
power, STt^ by ([Bond 1995]) 



exp(-wV2)W^] 



where 



^+1/2 
^o + l/2' 



(26) 



(27) 



10. The rms Amplitude Analysis 

It is sometimes convenient to characterize the fluc- 
tuations by their rms amplitude, A^^j,, given by 
^sfcy = ^totai - ^Lsf The raw rms of the data 



including instrument noise is Atotai = \^J2 ^T^/N, 
and the estimated contribution due to instrument 
noise is Ainst- 

The instrument noise, Ainst i cr^.^^^, is generated 
from the data covariance matrix, Cd, 4096 examples 
of data sets described by Cd are generated by 



(28) 



where N is a vector of the same dimension as t com- 
prised of samples drawn from a unit normal distribu- 
tion. The quantities Amst and cTAi„^t found from 
the distribution of the variances of the fake data sets. 
The uncertainty in Asky due to measurement noise is 
given by 

_mst Ainst 



^sky 



(29) 



\sky =c cr^;^:^ gives the am- 



This determination of A^^j, ± cr™** 
plitude (and uncertainty) of the fluctuations only over 
the region of sky which was measured, and is there- 
fore only an estimate of the amplitude of the overall 
distribution from which this sample was drawn. This 
is appropriate for comparing different measurements 
of the same sky but does not include sample variance. 
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The uncertainty in Asky due to sample variance 
((t£^^ ) is generated from a theory covariancc matrix, 
Ct, in the same way as cr/:^.^^^ was generated from Cjy; 
4096 fake data sets described by Ct, are generated by 



(30) 



The uncertainty cr^\ is found from the distribution 
of the variances of the fake data sets. This estimate of 
the uncertainty due to sample variance is larger than 
what is found by a Likelihood analysis; the latter in- 
cludes the effect of spatial phase information inher- 
ent in the model which the former does not (Knox 
IQSjll). The uncertainties and aX'*^* are added 

in quadrature to estimate the total uncertainty in 
Asky ■ The rms amplitude, Asky , and its uncertainties 
may be converted to STi using equation (16). 



11. 



Determination of 

The Frequency Spectral Index /3 



The first issue addressed in the analysis is whether 
the signal is due to the CMB or is due to some other 
foreground contaminant. As noted, the spectral in- 
dex of most foregrounds is expected to deviate from 
that of the CMB, f3 = 0. Since data are acquired 
in 6 frequency bins between 26 GHz and 46 GHz, /3 
for the fluctuations may be evaluated using the likeli- 
hood analysis described. Since the data acquired with 
the Ka-band radiometer only has adequate sensitiv- 
ity for the 3pt, 4pt, and 5pt synthesis, this may only 
be done internally for angular scales below £ « 100. 
The frequency span given by the Q-band data alone 
is inadequate to determine /?. 

To evaluate f3 the data from Ka93 3pt and SK94 
& SK95 3pt to 5pt are used. The covariance matrix, 
M in equation (||) has dimension^ 1422 x 1422. For 
the theoretical spectrum (eq. |2^) is used with m = 
0. The hkelihood analysis yields STi^ = 47^5 /xK 
and P — 0.21q3 for if, = 73. Repeating the analysis 



with m ~ — 1 yields /3 



0.2toi, and m 



-f 1 yields 



P = 0.1^0^3. Thus the limits on /3 are independent 
of choice of angular spectral index, m. The channel 
to channel calibration uncertainty is 2% to 5%. This 
causes an additional uncertainty in /3 of ±0.1, which 
is to be added in quadrature to the quoted errors. 



^ Ka93 has 3 frequency channels, 21 RA bins, and both East 
and West data (126 rows), Ka94, Q94 and Q95 each have 3 
frequency channels, 24 RA bins, and East and West data for 
3pt, 4pt and 5pt data (1296 rows). 



The spectral index, /3, of the fluctuations is con- 
sistent with that of the CMB {(3 = 0), and inconsis- 
tent with known potential foreground contaminants. 
In particular, contamination by Galactic synchrotron 
(/? = —2.8), Galactic free-free emission {(3 = —2.0) 
and Galactic dust (/3 = 1.7) in this data set are ruled 
out at large angular scales. Foregrounds are discussed 
in detail in Section |l^. 

12. Comparison With MSAMl-92 

Perhaps the most powerful systematic check avail- 
able with the Saskatoon experiment is the comparison 
with the MSAMl-92 experiment. MSAMl-92 is a bal- 
loon born experiment which observed the sky with a 
3pt beam at a declination of 82° in w 81 fields between 
14.4 hr and 20.4 hr RA ([Cheng ct al. 1994|). While 



the Saskatoon experiment observes the sky between 
26 GHz and 46 GHz, the MSAMl-92 experiment has 
observing bands sensitive to the CMB at 180 GHz and 
240 GHz, providing very wide frequency coverage be- 
tween the two experiments. The comparison of the 
two experiments of very different nature places strin- 
gent limits on systematic errors or foreground con- 
tamination. 

As mentioned in Section ^, data from 1995 which 
was acquired in the Ring mode (at a declination of 
82.05°) was multiplied by a weighting vector which 
closely synthesized the MSAM 3pt beam. To gener- 
ate this weighting vector, the elements Wi are adjusted 
to minimize the variance between the synthesized SK 
and MSAM beam patterns, assuming that the dec- 
lination of the two experiments is equal. In reality, 
the mean pointing of the two experiments differs by 
0.14°, and the MSAMl-92 experiment did not track 
in fixed elevation. This discrepancy is ignored in the 
in the analysis of A, but included in the Likelihood 
analysis. Once a SK weighting vector has been pro- 
duced which maximizes spatial overlap, it is scaled so 
that y^I{W) = 1.15, the value for the MSAM window 
function. To check the level of overlap between the 
the two experiments, the window function of the dif- 
ference of the synthesized beams is formed. For the 
difference beam y^I{W) = 0.25. Given this overlap, 
it is expected that over 90% of the sky signal will be in 
common between the two experiments. The window 
functions for the MSAM overlap are given in Figure |[ 

Given the high degree of overlap expected, a naive 
test is adequate to verify that the two experiments 
see the same signal. Figure ^ shows the two data 
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200 400 600 800 



Fig. 5.— The MSAMl-92 and SK95/MSAM Overlap 
window functions. The window functions are sUghtly 
different due to unequal beam size. The window func- 
tion for the difference is also presented. 

sets, their mean and difference. The reduced 
the sum is 3.43, while the reduced ^oi the differ- 
ence is 1.05, neglecting the anomalous third sample at 
14.58 hr in the MSAMl-92 data set. It is clear that 
both experiments see the same signal on the sky. 

To find the consistency of the amplitudes measured 
by the two experiments, an analysis of A as described 
in Section |l^ is performed. The data covariance ma- 
trix for the MSAMl-92 data {Cdms ) is taken to be di- 
agonal, neglecting correlations which were introduced 
by the removal of offset drifts in the data reduction. 
For the SK data, Cdsk made from the distribu- 
tion of the 18 second averages, and includes all known 
correlations. For analysis of the combined MSAMl- 
92/SK95 data, the weighted mean is used, and the off 
diagonal elements on the covariance matrix are gen- 
erated using 



SK 



•"MS 



){W 



(31) 



SK 



-"MSI 



where ' 



^SK 



and ' 



'MS 



are the statistical weights of the 



SK and MSAMl-92 data respectively. 

The results of this analysis are presented in Ta- 
ble ||. The errors on A include only cf"^^^* and do not 
include sample variance. This is appropriate for com- 
paring the consistency of two measurements of the 
same sky. The amplitudes are consistent, with SK 
favoring a higher value than MSAM at the la level. 

This analysis has not included calibration uncer- 
tainty in the two instruments. For Saskatoon it is 
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Fig. 6.— MSAMl-92 and SK95/MSAM Overlap com- 
pared. The first panel gives the MSAMl-92 3pt data. 
The second panel gives SK95/MSAM Overlap data 
for the same region. The third shows the mean and 
the fourth the difference. The data are in 85 equally 
spaced RA bins (only 81 have data), and are given in 
thermodynamic temperature units. The anomalous 
third bin in the MSAM data set is dropped in the 
analysis of the comparison. The 0.5° FWHM beam 
spans 0.24 hr of RA. 
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Table 8 
MSAMl-92 AND SK compared 



MSAMl-92 SK Mean Difference 

X^/v 2.66 1.82 3.43 1.05 

Atot 88.9 AiK 104.4 /zK 76.6 52.2 

Ainst ± crA,„,t 70.5 ± 5.3 79.5 ± 7.9 /iK 49.3 ± 3.9 /iK 53.2 ± 4.3 ^iK 

AsfeyicrX'/L 53.5 ± 7.0 /L(K 65.9 ± 10.2 /LtK 58.6±3.3^K 

STe{Asky) 46.5 //K 57.3 51 /iK 

(JTf (Full Likelihood) 52+^ VK 61^^? ^iK 50.9t6VK 



Note. — A comparison of the MSAMl-92 and SK95/MSAM Overlap data. The reduced entries are for the 
80 bins in Figure ^, neglecting the anomalous third sample from MSAMl-92. The uncertainties on A^^j^ do not 
include sample variance. The uncertainties on STi do. STi{Asky) is 5Ti inferred from A^^.j^ using equation (p^. 
Calibration uncertainties are not included. 



14%, and for MSAMl-92 it is 10%. Adding these in 
quadrature yields a relative calibration uncertainty 
of 17%. A least-squares fit to find the best rel- 
ative calibration between the two data sets finds 
Nsk/Nmsam = 0.82 ± 0.16. The relative calibra- 
tion based on the data favors a smaller normalization 
for SK or a larger normalization for MSAM, but is 
still consistent with the absolute calibrations used. 

The flat band power is found by a full likelihood 
analysis using the angular spectrum in equation ([2^). 
These results are presented in Table ^ The results 
for the two experiments are fully consistent. Note 
that calibration uncertainties have not been included 
in quoted errors. 

The original MSAMl-92 analysis is done in terms 
of the GACF, equation (|2^). This analysis is repeated 
here. For MSAMl-92, we find ^C{<d) = 70tg^/xK , 
using 9c = 0.92, which is consistent with the anal- 
ysis quoted in Chen g et al. (1994). (1994)[| For 
SK we find ^JC{0) = 81^1^ nK . Converted to 
band power estimates using equation (^6|) we find 
53 /iK and 61 jiK for MSAMl-92 and SK respectively. 

A determination of [3 of the fluctuations yields 



■iCheng et al. (1994) find 50 /iK < ^C(O) < 90AtK for the 5% 
to 95% confidence interval. Reanalyzing the MSAMl-92 data, 
we find 55 fiK < y/C{0) < 99 fiK . 



P = -0.1±0.2. The combination of the SK95 and 
MSAMl-92 experiments has essentially measured the 
frequency spectrum of the anisotropy of the fluctua- 
tions from 36 GHz to 240 GHz and have found them 
to have the frequency spectrum of the CMB. Over 
this large frequency range no other known source has 
this spectrum. Free- free emission has (3 — —1.45 (as 
compared to -2.1 in the Ray leigh- Jeans region) and 
interstellar dust has (3 — 2.25 (1.7 in the Rayleigh- 
Jeans region). The spectral index of these fluctua- 
tions is highly inconsistent with these foreground con- 
taminants. 

The differences between the two experiments ren- 
der it very unlikely that the observed signal is the 
result of mutually common contamination. MSAM is 
a balloon-borne experiment which observed the sky 
for several hours, while the Saskatoon experiment is 
ground based and observed the sky over many days. 
MSAM is a bolometer based experiment which ob- 
serves the sky at 180 GHz and 240 GHz while SK95 is 
based on 40 GHz HEMT amphflers. MSAM utihzes 
multi-mode optics while those of the Saskatoon tele- 
scope are single mode. MSAM calibrates on Jupiter 
while SK calibrates on Cas-A. The consistency of the 
measurements between MSAM and SK indicates that 
both experiments have measured the anisotropy of the 
CMB. 
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Fig. 7. — The angular spectrum in 37 bins. The re- 
sults for each synthesized beam for SK93-SK95 is pre- 
sented separately. These results are given in Table ^. 
The error bars do not include the ±14% (Icr) overall 
calibration uncertainty. 



13. Determination of 

The Angular Power Spectrum 

The angular spectrum is presented in three ways. 
The first is to find STi^ for the data set associated 
with each of the 37 synthesized beams, Hi , described 
in Section The results of this are presented in Ta- 
bic | and Figure |^. Also included in Table ^ are re- 
sults from an analysis of A as described in Section [Tc| . 
Note that for SK95 3pt through lOpt, the estimated 
contribution to the uncertainties due to sample vari- 
ance, crgXi' is approximately equal to the contribu- 
tion due to instrument noise, (Jst^- Consequently, to 
significantly reduce the error bars, both better sensi- 
tivity and more sky coverage will be required. This 
representation has significant spatial correlations be- 
tween entries. For example, the 3pt beams for Q95, 
Q94 and Ka94 all observe the same sky and at a sim- 
ilar angular scales. It is therefore inappropriate to 
combine entries simply by taking the weighted mean 
without taking into account the spatial correlations. 

The second way the data are presented is to com- 
bine data from spatially correlated synthesized beams 
between years. To do this, a covariance matrix 
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Fig. 8. — The angular spectrum grouped by spatial 
symmetry. These results are given in Table noL 
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Table 10: The angular spectrum grouped by spatial 
symmetry. The 14% overall calibration uncertainty is 
not included in the errors. 
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Table 9 

The angular spectrum in 37 bins 





Bins 


Atotal 


Ainst 






STe{Asky) 


—inst 
^5Tt 






6Te 






AQ iiK 
4y ^iv 




'?9 iiK 

OZ ^IV 


1 n4 

± . U4 


1 /xlv 


O ^iv 


1 1 riiv 


00-19 


'^fi + 12 iiK 

OU_g fll\. 


IC OA '^r.f 




DO 


AA 
44 


zLQ 
4y 


1.10 


49 
4Z 


A 
4 


1 9 
IZ 


00-18 


cn+lS 


4pt 


24 


48 


42 


24 


88 


97 


q 


w 


7P+21 
'0-19 


00+13 

OZ_g 


opt 


z^ 




oo 


9R 
zu 


7"^ 


OU 


1 n 

lU 


1 n 

lU 


96+21 
90-19 


44+17 


up h 


4:0 






oo 


fi9 
u. uz 


^4 

Ot: 


1 9 


Q 

y 


110-19 


oq + 16 


1 p L 


4:0 


SI 

1 


7^ 
/ o 


ou 


^4 


OU 


97 
z / 


lO 


1 '^4+22 


'^'^ <r 7'i 

oo , \ I o 


opL 


40 


111 

± ± ± 


lUl 


47 

4 / 


n 4R 


47 


ou 


1 9 
1 z 


1 C9+23 
10^-18 


Ou, ^ izy 


ypL 


zlS 
40 


1 97 


1 zlfi 


n 

u 


n 4^^ 


u 




1 9 
1 z 


l'^J-17 


fi ^ ^AP. 

U, \ 14U 


v^y^ opL 


Z4 


ou 


7"? 
i o 


4'^ 
40 


1 99 


47 


1 n 

lU 


1 1 
1 1 


oy_2i 


41-12 


4:pL 


9A 

Zi4 


oo 




OO 




^7 

O 1 


Q 


1 n 

1 u 


01+16 
°l-24 


19 ^ 7fi 

IZ, ^ ( u 


OpL 


9zL 

Z4 


oy 




99 
zz 


77 


9R 
zo 


1 Q 


1 n 

lU 


100+1^ 
1UU_23 


17 ^8 
1 ( , ^ 00 


UpL 


/IS 

40 


7^ 
1 o 


89 
o^ 


n 


fi7 



u 




Q 

y 


1ZU_24 


9^ <' fi9 
zo, \ UZ 


1 pL 


AR 

40 


1 OQ 




4^ 

40 




77 


97 

Z J 


1 9 

1 Zi 


ioy_24 


fi7 ^ 1 90 
U ( , \ 1 zu 


OpL 


40 


1 1 9 


1 07 


"^9 

OZi 


u.oo 


fil 
u ± 


ou 


1 1 
1 1 


1^7+19 

23 


fi4 ^ 1 1 q 

U4, ^ iiy 


ypL 


AR 

40 




1 9Q 


44 

44 


U.'lO 


Ql 


oo 


1 1 

1 1 


1 '0-23 


-,00+59 
10S_43 


v^yo oci// opL 


9zL 

Z4 


1 01 


78 
/ o 


DO 


1 9Q 


4Q 


c 
o 


1 1 
1 1 


64+12 


00-14 


^pL 


9zL 

Z4 


DO 


ou 


OO 


1 09 

_L .UZ 


"^9 

OZ 


1 1 

1 L 


1 n 

lU 


00 + 13 
°'J-29 


00 + I6 
0»_8 


OpL 


9A 
Z4 


fi9 


40 


44 
44 


n SI 
u.oi 


oo 


u 


Q 

y 


1 ns+8 

iU»_32 


C4+I9 


UpL 


40 


7fi 


fi8 

OO 


4Q 
4y 


n 71^ 


1^4 
u^ 


e 
o 


s 
o 


1 '^f^+^ 

100_37 


fi4+18 
04-11 




4R 


74 




^9 


67 


1 o 


7 


1 


100_3g 


79 + 19 
'^-12 


opt 




7f^ 


oo 




u.Oc/ 


R4 
o^ 


Q 

y 


1 n 

lU 


1 78+^ 


*1-12 






Uc/ 


"^4 


4^^ 


^4 


7Q 


1 
1 u 


Q 

y 


ly ' -37 


76+2" 

'0-13 


lOpt 


48 


70 


60 


37 


0.50 


74 


16 


11 


217+L 

' -38 


76+?^ 

' "-16 


Upt 


48 


63 


56 


31 


0.46 


66 


18 


11 


237t«7 


83t?2 


12pt 


48 


80 


60 


53 


0.42 


126 


13 


11 


9^7+8 
Z0(_3g 




13pt 


48 


78 


59 


51 


0.40 


127 


14 


11 


277+1^ 




14pt 


48 


74 


64 


38 


0.37 


102 


24 


11 


297+1^ 


74+26 
'4-24 


15pt 


48 


78 


73 


28 


0.35 


79 


48 


10 


316^Ji 


1 94+35 
1Z1_26 


16pt 


48 


96 


73 


63 


0.33 


190 


22 


8 


333^^!^ 


1 1 0+38 
110-38 


17pt 


48 


95 


78 


53 


0.31 


172 


30 


8 




110, < 197 


18pt 


48 


97 


84 


48 


0.29 


166 


41 


8 


ci8/!_3Q 


88, < 191 


19pt 


48 


99 


91 


40 


0.27 


150 


63 


8 


404+13 


0, < 163 


Q95 Ring 3pt 


96 


107 


77 


74 


1.14 


65 


5 


7 


170+^^ 


59+11 


4pt 


96 


104 


78 


68 


0.85 


80 


8 


9 


234^67 


mi 


5pt 


96 


101 


89 


48 


0.70 


69 


18 


8 


286+39 




6pt 


96 


106 


107 





0.59 







8 


332^99 


19, < 90 



Note. — Results from each synthesized beam from SK93, SK94 and SK95. The SK93 entries are from Netterfield et al. 
(1995). The Ka94 and Q94 entries differ slightly from those found in Netterfield (1994) due to the inclusion of East to West 
correlations in the analysis. l\totai is the raw rms of the binned data. l\i„st is the rms predicted by the data covariance 
matrix, Cd- A^sky = A^otai ~ ^nst gives the sky rms. 5Ti{/^sky) gives the band power inferred from /S.sky using equation 
(0)- ""iTf* estimates the contribution to the error bars due to the instrument noise and of^^ estimates the contribution to 
the error bar due to sample variance. The limits on are for where the window function falls to e~^^'^ of the peak value. 
The final column lists JTf from the full likelihood analysis. The overall 14% caUbration uncertainty is not included in the 
errors. 
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Table 11: The angular spectrum in 5 bins. The first 
bin is comprised of the 3pt to 6pt Cap data from all 
three years. The second is from the 7pt to 9pt Cap 
data and the 3pt Ring data. The third is lOpt to 12pt 
Cap and 4pt Ring . The fourth is 13pt to 15pt Cap 
and 5pt Ring . The fifth is 16pt to 19pt Cap and 6pt 
Ring . See Section ^ for a description of the cut levels, 
m is the angular spectral index within each bin. 



Cd + Ct is made which includes all of the data corre- 
sponding to a given spatial symmetry. For example, 
the Ka94, Q94 and Q95 3pt data forms a 144 x 144 
matrix. The band power is then evaluated using a 
maximum likelihood test. An effective window func- 
tion is generated using equation ([2^). The Ka93 3pt 
data has largest spatial overlap with the 94/95 4pt 
data and is included there. These results are pre- 
sented in Table |l0| and Figure ^. The tendency for a 
rising spectrum is evident here though the scatter in 
the points, particularly at high is still considerable. 

Finally, the data are combined into 5 angular scale 
bins and presented in Table 



11 



and Figure To pro- 
duce each of these bins, all of the data which is sen- 
sitive to a given range of angular scales is combined 
using the full correlation matrix (see Section 0) . For 
instance, the first bin is comprised of the 3pt to 6pt 
Cap data from all three years. The effective window 
functions for each bin using equation ( |2^ ) are pre- 
sented in Figure |l^. 

Also presented in Table ^ are values for the an- 
gular spectral index m (see Section The negative 
spectral index around £ « 80 which was reported in 
Netterfield (1994) with the 1994 data alone is again 
seen in the lowest angular bin. The lowest four points 
in Figure || also shows this locally falling spectrum. 
However, this trend is dominated by the 3pt data be- 
ing high. Additionally, for each of the 5 angular bins, 
including the first one, the slope is consistent with 0. 
The angular spectrum of the whole sky may have a 
dip at these scales, or this effect could be a result of 
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Fig. 9. — The angular spectrum in 5 bins. The results 
are tabulated in Table |ll| and the effective window 
functions are shown in Figure |l^. Predicted spectra 
from six representative theories are also shown. They 
are, from top to bottom &i I — 160 a flat A-I-CDM 
Model with SIa — 0.7 (Ratra and Sugiyama 1995), 
Standard CDM (Ratra et al. 1995), a Texture model 
(Crittenden and Turok, 1995) a PPI model (Peebles, 
1995) an = 0.4 open model (Ratra et al. 1995), and 
a model with reionization (Sugiyama 1995). The tex- 
ture model is arbitrarily normalized. All others are 
normalized to COBE. To accurately compare each of 
these spectra with the data, they must be convolved 
with the window functions in Figure ^ The 14% 
overall calibration uncertainty is not included in the 
error bars. This effects the normalization of the spec- 
trum, but not the shape. 
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Fig. 10. — The 5 bin effective window functions. 
These have been generated from the window func- 
tions of the contributing effective antenna patterns 
using equation (122" 



sample variance. This issue will be resolved by ex- 
tended sky coverage to reduce sample variance. 

While all known correlations have been taken into 
account in generating the likelihoods for each of the 
combined groups, there are still some correlations be- 
tween groups. They are small (< 0.2), however, and 
may be ignored for most analysis. 

Not included in the quoted error bars are the effects 
of calibration and beam uncertainties. As discussed, 
the calibration uncertainty is 14% (Icr), mainly due to 
uncertainty in the temperature scale of Cas-A. This 
has very little effect on the uncertainty of the spec- 
tral index, /?, and no effect on the shape of the angu- 
lar spectrum. Rather this has only an effect on the 
over-all temperature scale. This large normalization 
uncertainty must be included in any comparison with 
theoretical predictions. 

Beam uncertainties, as well as contributing to cal- 
ibration errors, also produce uncertainties in the an- 
gular spectrum. For uncertainties in the beamwidth 
along the elevation axis, the only effect is in the cali- 
bration (miss-estimating the beam solid angle in pre- 
dicting the signal from Cas-A) , and has been included 
in that error bar. For uncertainties in the azimuth 
beam width the effect is a function of the weight- 



ing vectors used. For the 4pt synthesis, where the 
physical beam width is small compared to the width 
of the synthesized lobes, only the temperature scale 
is effected, as with the elevation beam. For the 13pt 
synthesis, however, the additional effect of mis-stating 
the beam overlap between points very nearly cancels 
the effect on the calibration. The result of this is 
that the 2% azimuth beam uncertainty yields an ad- 
ditional 2% relative uncertainty between the 4pt and 
13pt synthesized beam data. This effect is small for 
this experiment given the error bars on the spectrum 
and is ignored. 



14. Foreground Contamination 

The angular spectra presented in the previous sec- 
tion are for the microwave sky at 40 GHz. There are, 
however, potential sources of foreground contamina- 
tion. Two major classes are diffuse Galactic emission, 
and unresolved point sources. Limits may be placed 
on the level of both of these potential contaminants. 

Potential sources of diffuse Galactic contamination 
include dust, synchrotron, and free-free emission. The 
level of dust contamination at 40 GHz and £ = 100 
has be estimated to be less than 2/iK (Tegmark & 



Efstathiou 1995), which is far smaller than the sig- 
nals measured here. Similarly, recent measurements 
of the amplitude of diffuse free-free emission based on 
Hq emission have been made for the Saskatoon ob- 



servi ng region ( Gaustad et al. 1995 ; Simonetti et al 
1995| ). These measurements place an upper limit of 



6 fiK at 30 GHz for the SK93 3pt beam. Once again 
this is too small to account for the observed signals. 
An analysis of high latitude Galactic emission in the 
COBE 2 year data places similar limits on both fore - 
grounds at large angular scales ( Kogut ct al. 1995| ). 
An estimate of the amplitude of synchrotron emission 
can be made by extrapolating radio maps at 408 MHz 
(iHaslam, C.G.T. et al. i982| ) and 1.4 GHz ( [Reich & 



Reich 1986) to 30 GHz assuming a temperature spec- 
tral index of /3 = —2.8. This predicts a signal of 
less than 5 /iK in Ka93, which is much smaller than 
the measured signals. However, if the spectral in- 
dex of the synchrotron emission varies spatially by 
SP « 0.05, then 30 /zK signals could be expected. 
The spectral index of these fluctuations would still be 
characterized by /3 « —2.8. 

The most powerful discriminant against foreground 
contamination is based on the frequency spectrum of 
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Table 12 
Bright Point Sources 



Source 




Flux at Vo 


Spectral Index a 


Flux at 40 GHz 


Reference 


0014+81 


10.7 GHz 


0.73 Jy 


0.36 


1.2Jy 


S5 


0210+860 (3C61.1) 


30.0 


1.0 


0.2 


1.1 


HR 


0454+844 


10.0 


1.5 


0.0 


1.5 


K 


0615+820 


10.7 


0.86 


0.0 


0.86 


E,K 


0740+82 


10.7 


0.65 


-0.47 


0.35 


HR, S5 


1003+83 


10.7 


0.66 


-0.10 


0.58 


S5 


1039+811 


90.0 


0.80 


0.0 


0.80 


IRAMc 


1053+815 


20.0 


0.4 


0.43 


0.54 


E,VLBIc 


1050+812 


10.7 


1.1 


-0.1 


0.96 


K 


1221+809 


10.7 


0.6 


+0.2 


0.78 


S5, VLAc 


1637+8239 (NGC6251) 


10.7 


0.8 


-0.4 


0.47 


HR 


1637+826 


10.7 


0.73 


-0.4 


0.43 


S5,VLAc 


2342+821 


30.0 


0.8 


0.2 


0.85 


HR, S5 



Note. — Listed are parameters for the brightest point sources used in the estimate of the level of point source 
contamination. The flux at 40 GHz is extrapolated from the measurement at by S^v) = S'(fo)(^'/^'o)", where 
a = /? + 2. The measured value from the literature for a used where possible. Otherwise a = (flat spectrum or 



constant flux with frequency) is assumed. The data are obtained from Eckart et al. 1986 (E), Herbig & Readhead 
19^ (HR), ^teppe ct al.T992| (IRAMc), |Kiihr et al. 1981a] (K), |Kiihr et al. 1981b| (S5), [Perley 1982| (VLAc) and 
Edelson 1987| (E). 
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Fig. 11. — Estimated levels of foreground contami- 
nation at 40 GHz. The 5 bin data from this experi- 
ment are given as asterisks. The estimated angular 
spectra for the various foreground sources are from 
Tegmark and Efstathiou (1995). The free-free angu- 
lar spectrum has been renormalized for our observing 
region using a limit based on Ha emission (Gaustad 
et al. 1995; Simonetti et al. 1995). An upper limit 
for the normalization of the synchrotron spectrum has 
been estimated by extrapolating the 408 MHz Haslam 
map to 40 GHz for our observing region (Haslam et 
al. 1982, WoUack et al. 1993). The square boxes rep- 
resent contributions to our data estimated by extrap- 
olating exiting point source surveys to 40 GHz (see 
text). Also indicated is an upper limit on sources 
at smaller scales based on the OVRO Ring experi- 
ment (Meyers, Readhead, and Lawrence, 1993) ex- 
trapolated from 20 GHz. 



the fluctuations (/?). In Section 111], for ^ « 73, /3 
was found to be 0.2 ± 0.3 which is inconsistent with 
the spectral index of dust {J5 = 1.7), free-free emis- 
sion [jj = —2.1), or synchrotron emission (/3 = —2.8). 
Additionally, since the angular spectrum for galactic 
emission falls as 6Te cx ^"^/^ ( fTegmark fc Efstathiorj 
"T995| ), it is very unlikely that the galaxy is a signifi- 



cant contaminant at any angular scale. 

Of more concern for several reasons is potential 
contamination by unresolved point sources. Firstly, 
the angular spectrum of a family of point sources is 
expected to rise as STi cx £, so the f3 limits placed 
at large angular scales are weaker at smaller angular 
scales. Additionally, since the mechanism for radia- 
tion varies between individual sources, the expected 
(3 also varies, and in many cases is not accurately de- 
scribed by a simple power law. This is of particular 
concern since no flux-limited survey of the SK observ- 
ing region in has been performed over our frequencies 
(26 GHz to 40 GHz). However, it is still possible to 
make some estimates of the level of contamination 
via comparison with MSAMl-92, source counts from 
the OVRO Rin g experiment ( Meyers, Readhead, & 



Lawrence 1995 ) at smaller scales and an extrapola- 
tion of existing point source surveys. 

While (3 varies from source to source, it is expected 
that over a sufhciently wide frequency range no fore- 
ground source will have a spectrum consistent with 
(3 = 0. Consequently our comparison of SK95 with 
MSAMl-92 which yields a hmit on /3 of 0.1 ± 0.2 with 
a frequency baseline from 40 GHz to 240 GHz assures 
that the sky covered by both experiments is devoid of 
significant point contamination. However, this field 
only represents 8 square degrees of the total of 200 
square degrees observed by SK95. 

A second approach is to use results from other ex- 
periments and point source surveys. The OVRO Ring 
experiment has observed 96 fields at 20 GHz. They 
find STi = lOi^iK at £e = 2017. Here we make the 
assumption that the signal is dominated by Poisson 
distributed unresolved point sources. Extrapolating 
in angular scale by STi oc £ to £ = 404, the angular 
scale of the 19pt data set, one arrives at an expected 
contribution of 21 . Since OVRO Ring was per- 
formed at 20 GHz, this number needs to also be scaled 
in frequency. A conservative estimate of the 'typical' 
(3 for a point source is that its flux is constant in 
frequency, (i.e., flat spectrum or /3 = —2). Thus, an 
estimate of the point source contribution at the small- 
est scales of the SK experiment based on OVRO Ring 
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source counts is 10 /iK at 40 GHz. 

As a final estimate of tlie level of point source con- 
tamination, the expected amplitude of known point 
sources in the observing region is extrapolated. Ta- 
ble |l^ lists the brightest of these. In addition to the 
sources list ed in Table [l^, all of the sources in the 
S5 catalog ( Kiihr et al. 19811: ) are used. Based on 
an analysis of IRAS faint sources, it is considered un- 
likely that a new population of sources exists within 
our observing frequencies (Condon et al. 1995). The 
spectral index from the literature is used to extrap- 
olate to 40 GHz whenever possible. If a measured 
frequency spectral index is not known, the source is 
assumed to be of constant flux with frequency. A 
source map of the observing area is made, and the 
SK observing strategy simulated on it. The largest 
estimated contribution is to the 19pt data, and is 
5Ti = 35 /iK , which could be significant. However, 
from Table ^, the 19pt data yields a 163 /xK upper 
limit with a maximum likelihood at fiK . The esti- 
mated contribution to the 17pt data is 23 /iK , which 
would imply a 4% contribution to our 69 /iK detection 
at these scales. The predicted signal is dominated by 
sources in RA bins 5 and 10. The measured 17pt 
data show no excess in these bins, and re-analyzing 
the data with these bins deleted changes the most 
likely value for the 17pt data by less than 1%. For 
the 13pt data, a,t £ — 277, the expected contribution 
is STi = 12 /iK . From Tabic ^, the measured value 
is 115 /iK . Subtracting the source contribution in 
quadrature yields less than a 1% effect. 

The data, along with estimated contributions from 
foreground contaminants are presented in Figure |l^. 
The indications are that the SK data are not seri- 
ously contaminated by foreground sources. However, 
extending the coverage to smaller angular scales will 
require improved discrimination from point sources. 

15. Discussion 

The rich set of internal and external consistency 
checks inherent to this experiment indicate that a sig- 
nal fixed to the celestial sphere has been detected. 
The limits placed on the frequency spectral index /? 
rule out significant contamination from diffuse Galac- 
tic emission. The limits on f3 over a very large fre- 
quency range placed by the successful comparison 
with the MSAMl-92 experiment further reduce the 
likelihood of significant contribution to the signal by 
foreground contaminants. Conservative extrapola- 



tions from existing point source catalogs strengthen 
this argument. The evidence is compelling that the 
signal which has been detected is due to anisotropics 
in the CMB. 

This experiment has measured the angular spec- 
trum of these anisotropics simultaneously with one 
experiment and with equal calibration error at all an- 
gular scales. We find that the spectrum rises from 
STi = 49^5 /iK at £ = 87 to 85^3" /iK at £ ^ 237. 
This can be used to place significant limits on the 
reionization history of the universe, and, within a 
given theoretical framework, begin to place limits on 
cosmological parameters. A quantitative comparison 
with theory will presented elsewhere. 

Given the relative complexity of this analysis, it is 
important to mention that much of the analysis has 
been repeated with completely different programs by 
two of the authors to verify accuracy. The consis- 
tency of the naive A based tests with the full like- 
lihood analysis further strengthens our confidence in 
our programs and techniques. In addition, Andrew 
Jaffe and J.R. Bond at CITA have repeated and ex- 
tended our analysis of the SK93 and SK94 data, again 
confirming our analysis. 

To improve the discriminating power of this exper- 
iment several things are needed besides reducing the 
size of the error bars. First of all, it must be confirmed 
at high i. Additionally, a wider frequency coverage 
at smaller scales will improve discrimination against 
point sources. A better calibration will improve com- 
parison with COBE, effectively increasing the range 
of angular scales over which a theory can be tested. 
It is possible to make a map of the CMB anisotropy 
from these data. A map of an extended region of the 
sky with better signal to noise will be useful to distin- 
guish between Gaussian and non-Gaussian models. 

We are indebted to Dave Wilkinson for his work 
on the experiment, and for many very useful discus- 
sions. We thank Marion Pospieszalski and Mike Balis- 
ter at NRAO for providing the HEMT amplifiers upon 
which this experiment was based. George Sofko and 
Mike McKibben at the University of Saskatchewan 
and Larry Snodgrass at the Saskatchewan Research 
Council provided valuable site support in Saskatoon. 
This work has benefited from useful discussions with 
Dick Bond, Ruth Daly, Tom Hcrbig, Andrew Jaffe, 
Jim Peebles, Bharat Ratra, Dave Spergel, and Paul 
Steinhardt. Financial support has been provided by 
NSF grant PH 89-21378, NASA Grants NAGW-2801 
and NAGW-1482, a Research Corporation Award, a 
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Packard Fellowship, and an NSF NYI grant to L. 
Page. 

The data, beam profiles, window functions and 
covariance matrices will be made publicly available 
upon acceptance of this paper. 
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